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Abstract 

In this work, we consider the London limit of the (2 + 1)D Georgi- 
Glashow model with dynamical quarks. Following Polyakov's monopole 
plasma approach and using dual methods to treat the relevant matter 
sector, we derive in a clear and straightforward manner the recently 
proposed gauging of the discrete Z(N) symmetry of the associated ef- 
fective vortex theory. Our procedure applies to bosonic as well as to 
fermionic matter, enabling to derive a useful representation for the 
Wilson loop and discuss string breaking, due to the creation of dy- 
namical quark-antiquark pairs. This phenomenon corresponds to a 
perimeter law, a behavior that has been already observed in the lat- 
tice. 



One of the most interesting problems we know at present in Physics corre- 
sponds to look for possible phases and understanding confinement in strongly 
interacting systems like QCD. Then, the search for alternative methods en- 
abling to describe at least the relevant variables that determine the possible 
phases in the nonperturbative infrared regime becomes a central problem. 
One of these alternatives corresponds to duality; fundamental progress in 
this direction has been developed many years ago by A. M. Polyakov and G. 
't Hooft pig. 
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In particular, the dual superconductivity scenario 'T\ can be realized in 
the framework of the pure SU(N) Georgi-Glashow model in (2 + 1) dimen- 
sions. Classically, this model contains vortices with topological charge Z(N). 
At the quantum level, the vortex sector can be represented by means of 
vortex operators associated with instanton (monopole-like) singularities in 
euclidean spacetime, where the vortices are created or destroyed. There are 
two relevant Green's functions: two-point and iV-point correlators. They are 
incorporated by means of the effective lagrangian for the vortex field j2], 

d^Vd^'V - fi^VV -\{VVf - P{V^ + V^), (1) 

which displays a global Z(N) symmetry (for a review, with further advances, 
see jS] and references therein). In the Higgs phase, when the vortex is an 
elementary excitation, the mass parameter /x^ is positive and there is no 
spontaneous symmetry breaking of the Z(N) group. Then, the vacuum ex- 
pectation value of the vortex field vanishes. However, when the vortices 
condense, the VEV of the vortex field is nonzero and Z(N) is spontaneously 
broken. This leads to domain walls between different possible vacua which 
for planar systems are string-like. When considering heavy probes, corre- 
sponding to a "quark-antiquark" pair, a string is formed in between thus 
leading to an energy which is linearly increasing with the pair separation, 
and an area law for the Wilson loop. This phase is then confining. 

This nice mechanism continues to be an important topic of research; it 
is a defying problem with regard to its possible implementation in QCD [2], 
|3], while many open questions are still present in simplified models, which 
are also interesting in a condensed matter context [3 E] ■ 

Confining systems in (2 + 1)D that incorporate matter in the fundamen- 
tal representation are not so understood. However, the Wilson loop still 
provides information about the vacuum response to the presence of a heavy 
quark-antiquark probe. In particular, a perimeter law for the Wilson loop 
would imply a constant potential, coming from the creation of dynamical 
light quark-antiquark pairs out of the vacuum. This phenomenon is called 
string breaking and has been observed in lattice QCD j7j (using the Polyakov 
loop) and lattice (2 + 1)D models containing bosonic or fermionic dy- 
namical matter. In particular, in ref. ^9j, string breaking has been detected 
by means of the Wilson loop criterium. 

In ref. ^U], by using the vortex condensation approach to the (2 + 1)D 
Georgi-Glashow model, the authors proposed the gauging of the Z(N) sym- 
metry in eq. ([[J, when bosonic matter in the fundamental representation 
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is included. They also argued that the perimeter law is associated with a 
boundary effect coming from this gauging. 

Besides the Georgi-Glashow model, confinement can be also understood 
in the case of pure compact QED(3) |jy. In fact, both models are closely 
related as, in the London limit, the relevant low energy dynamics of the 
former is an abelian subgroup of SU(N), which is compact. 

The aim of this letter is following the monopole plasma approach in- 
troduced by Polyakov to provide a simple derivation of the Z(2) gauging 
phenomenon when dynamical fermionic quarks are present. Also, upon ob- 
taining a dual representation for the Wilson loop, we will discuss how the 
perimeter law is attained in the limit where the quarks become massless. 

Then, our discussion will be focused on compact QED(3), eventually 
considered as the London limit of the SU(2) Georgi-Glashow model, when 
the Higgs and become very massive. Let us take the following correlation 
function. 



where is an external source to probe the system and = Efj^^ypd'^A^ is the 
dual tensor for an abelian gauge field A'^. The euclidean action associated 
with our model is. 



In the context of the Georgi-Glashow model, J'^ and refer to the SU(2) 
current and gauge fields along the U(l) internal direction which is unbro- 
ken before compactification, while ip transforms in the fundamental repre- 
sentation of SU(2) and is formed by a pair of four component Dirac fields, 
describing a parity symmetric system. 

In particular, note that the correlator (j2)) will be useful to study the 
Wilson loop, W = {expig § dx^^A^), integrated over a boundary 9S, as this 
loop can be rewritten, by using Stokes theorem, in terms of a (S) surface 
integral of the dual tensor. That is, W = W{g6j]n), where 5s is a Dirac 
5-distribution with support on S and is the field of unit vectors normal 
to S, defined by the surface element dS^ = nfj,dS. This is one of the steps 
we will need in order to compute the Wilson loop when compactifying the 
^-sector of our model, as this process is accomplished by first performing 







(3) 
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the replacement, 

r^r+f , nx)=gm [ dy^6^'\x-y), ggm = 27T, (4) 

where 7 is a Dirac string associated with instanton (monopole-hke) singulari- 
ties located at the string endpoints, and then summing over all the instanton 
anti- instanton configurations (monopole plasma). In this manner, processes 
where ^-vortices are created or annhilated are incorporated. The vortex 
magnetic flux is given by the monopole magnetic charge Qm, which has been 
taken to saturate Dirac's quantization condition. 

Then, we are left with the problem of implementing in the second term 
of the action Q. For this aim, it is convenient to use a dual transformation 
of the type introduced in refs. [TTl |T2I, and extensively discussed in refs. 

Kp + ig J d^xA^Jf, ^KB + igJ d^x A'^e^.^pd'^AP, (5) 

(see ref. ^3]) where the U(l) fermion matter current is bosonized using a 
vector fleld A^ whose gauge invariant action is given by. 

The functional T\6\ = — lndet(^ + m + p) is the effective fermionic action in 
the presence of an abelian gauge fleld 6^, along the abovementioned unbro- 
ken U(l) direction. In this regard, note that if we exponentiate the second 
member in (jS)) and then use eq. ©, we can path-integrate over A to ob- 
tain a 5-functional leading to the constraint, 6^ = gA^ + pure gauge. Then, 
integrating over we are left with exp(— r[5f^]), that is, the fermion path- 
integral using the flrst member in (jSJ. Then, the mapping deflned by eqs. 
(0) and ^ is exact. The current bosonization rule <-> e^ypd'' A^ is also 
universal, and it has been used to discuss the universal transport properties 
in quantum Hall systems (see ref. P3])- This mapping generalizes the well 
known bosonization of (1 -|- 1)D fermionic systems, as in that case (0) and 
(jU)) are valid with the replacement epupd^A^ — >• Epud^cp, that is, the usual 
bosonization rule for the fermion current, see also ref. [T3] . 

Now, in the (2 + 1)D case, because of the current bosonization rule, we 
can use an integration by parts in the second member of (0) to rewrite the 
interaction in terms of f^, thus enabling a proper implementation of (jlj). 
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Then, in compact QED(3) with dynamical fermions, the correlator Q can 
be represented as, 

W{s) = VjVAVAe-^+'f'''^'>^'^^'+^'\ (7) 

S = jd^x]^{f+jf + igjd'x{f + j)A + KB[A]. (8) 

The measure Vj represents the path-integral over the monopole plasma, 
whenever the integrand only depends on the monopole density. See ref. 0, 
for another possibility to implement the compactification in QED(3) models 
with dynamical matter; implications associated with both possibilities are 
interesting in the context of applications to condensed matter systems E] ■ 
In our case, we can follow the analysis for pure compact QED(3) given 
in ref. ^Tj, adapted for the computation of the correlation function W{s). 
Introducing a Lagrange multiplier to linearize the Maxwell term in eq. (jH)), 
and integrating over the field A^, we obtain the constraint, gA^ — — = 
d^cp. Using this information in eq. (jS)), we get, 

W{s) = ^J l?jPAD0e-/^'"(^(^<^+^-^^)'+'^'^^'"^)-^^[^l, (9) 

where we integrate over 0, instead of A. Then, the second term can be 
integrated by parts and written in terms of the monopole density p = d^j^^, 
so that the integration over j is done over the monopole plasma, thus leading 
to a sine-Gordon term, cos gmfp DO- After the field redefinition gm(f) ^ (f), 
and using Dirac's quantization condition ggm = 27r, we get, 

1 ~(d?x(-^{d<t>+gn^s-2TTAf~ifi^cos<p\~KB\A] 



W{s) = J VAV(f)e ' V^.^^"'"'" ' \ (10) 

where ^ is the monopole fugacity and /i is a constant with the dimension of 
mass. Of course, if matter were not coupled, eq. (fTO|) would read, 

Wo{s) = ^ jvAV^Pe-^''^^'''^'-''^-'''-'). (11) 

It is worth underlining here that Polyakov's dual field can be identified 
with twice the phase of the vortex field V (see j^] and references therein). 
Then, the symmetry (j) + 27rm in eq. (fTT|) corresponds to the Z(2) global 
symmetry of the N = 2 effective vortex theory in eq. ([Q). 
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When = gSY;n^, the correlator VFo("5) gives the Wilson loop for pure 
compact QED(3), here represented by means of the celebrated dual action for 
the field (see ref. ^J). According to Polyakov's analysis, the sine-Gordon 
term introduces a gap, leading to confinement of heavy probes. This comes 
about from a saddle point evaluation of eq. (fTTj) . 

df,{d>'^ + 27iS^n'') - a^sm^ = , = c/^^/il (12) 

Outside S, (p satisfies the sine-Gordon equation and is given by a domain wall 
containing a 27r discontinuity at S, coming from the source term. Therefore, 
when computing (dcj) + gms)"^ in eq. (fTTj) . the associated delta singularity 
in d(j) cancels against the = 27T6j:,n term. Away from S, the field is 
exponentially suppressed. Then, for a large surface, the domain wall leads 
to an almost uniform action density localized on E implying the area law for 
the Wilson loop, that is, confinement. 

Now, we see that in the case where the dynamical quarks are coupled, 
the theory is strongly affected. Namely, the global Z(2) present in eq. (Illj) 
becomes local. That is, the symmetry in eq. (fTUI) turns out to be, —>■ 
(p + 2nm{x), Afj_ Afj_ + dfj,m{x), where m{x) is an integer valued function. 

According to the discussion in ref. ^H]; the main operational difference 
between the local and global versions of a discrete symmetry amounts to a 
different prescription to carry out functional integrals. In the local case, the 
fields (f) and take values in the quotient space. Then, with the equivalence 
= + 271, must be considered as an angle. As the dual action i^sfA] 
is gauge invariant for any associated transformation A ^ A — (l/27r)50, 
we can decouple the field in the correlator (fTUI) . This also applies to 
the longitudinal component dip of the external source QmS appearing in the 
correlation function (jTU)), 

r d' s^{x') 
J A7r\x — x'\ 

whenever ijj corresponds to an "angle" field configuration. For the Wilson 
loop we have, gmd^s^ = 27in^d^6-£ and the values -0^, when we tend to a 
point on S from above or below, verify ip^ = ip~ + 27i, which are identified. 
Then, we can absorbe the dip term in W{s) by means of another transfor- 
mation of the field A, and the correlator (fTUI) becomes, 

W{s) = j^jVAe ^ ^ . (14) 
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As the system we are considering is parity symmetric (see below eq. 0), 
-^'^[yl] does not contain a Chern-Simons term; it only depends on the dual 
tensor F = edA. Then, transforming A —>■ A + R and using eq. (fT^ . the 
exponent in ()14p will contain the term / (Px tt^A? ^ plus the bosonized action 
evaluated on F + edR = F + Qm^ds (cf. eq. (O)- 

This leads to a qualitatively different situation when compared with the 
pure gauge result in eq. (fTT|) . There, the relevant source for the Wilson 
loop is concentrated on a surface, and the localization scale in the dual sine- 
Gordon theory implies the area law. On the other hand, when dynamical 
quarks are present, the local Z(2) symmetry changes the relevant source 
to be eds\^ = e'^'^Pupd^Sj:, which is tangent to the boundary of S and is 
concentrated there. For instance, in the case of a planar surface on the 
X3 = plane, the nonzero components are, +5'5(x3)920e and — (7 5(a;3)(9i9s, 
for /i = 1,2, respectively, where Gs = 1, when the transverse coordinates 
(xi,X2) are on S, and zero otherwise. 

This change in the dimensionality of the source indicates that, depending 
on the localization properties in eq. ()14|) . a perimeter law could be observed 
in the phase where dynamical quarks are coupled. 

In order to go further, let us analyze the general structure of the fermionic 
effective action in eq. (jH))- Again, because of parity symmetry, T[b] depends 
on the dual tensor G = edb, 

r[b] = ^ J d^xGO{-d^)G + nonquadratic part. (15) 

where the nonlocal operator 0{—d^) in the exact quadratic part comes from 
the exact vacuum polarization tensor. For large masses, it can be expanded 
according to, O = (l/m)(ci + C2 (^) + 03(^)2 + ...). The q's are adi- 
mensional constants with alternating signs. Because of charge conjugation 
symmetry, only terms with an even number of external legs contribute to 
T[b]. So that the nonquadratic part starts with four G factors whose leading 
term is 1/m^. In general, a finite order expansion of 0{—d^) presents un- 
physical poles, which are avoided by the implicit cut-off. This can be already 
seen at order 1/m^, as ci > and C2 < 0. However, this can be improved by 
keeping the full momentum dependence of the polarization tensor (see ref. 

m)- 

We will consider two important cases where additional information can 
be given: i) large m, where the nonquadratic part in F is suppressed, ii) 
m — > 0; in this case, using renormalization group arguments, it has been 
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argued that the infrared physics is associated with a nontrivial stable fixed 
point, dominated by the quadratic term in eq. (fT^ . where O ~ a (—9^)"^/^ 
(see ref. and refs. therein). For a related discussion see also ref. jl9j . 

Then, in these limiting situations, the (fixed point) quadratic fermionic 
effective action, together with eq. (jH)), leads to a quadratic bosonized action, 

Kb[A] = i I d'xF[-d^O{-d')]-'F, (16) 

and path integrating in eq. ()14|) . we arrive to an expression that can be used 
to study the infrared Wilson loop behavior. 



In case i), using O ~ (ci/m) ^ in eq. (fTTj) . we see that the exponent 
takes the form {l/2g'^)R'^. This result is physically sensible: for any finite 
probe separation, if we start at a certain mass scale where string breaking 
occurs, as the fermion mass increases, a situation is expected where the 
dynamical quarks cannot produce a complete screening anymore. This is 
manifested through the "magnetic" energy density form of the exponent, 
and the infinite localization scale in the field R, 

R = ed[-d'^]-^J , J = gmeds, (18) 

which prevents the onset of a perimeter law for any given finite loop. 

When m — i> 0, as discussed, we can still use eq. (fT^ to analyze in- 
frared behavior. For this aim, let us consider a Wilson loop, running along 
a large line of size L contained on the Xi-axis, which then crosses a trans- 
verse distance of the order of L, goes along a line parallel to xi and crosses 
back to close the loop. The part on the xi-axis gives a nonzero contribution 
eds\i = g6{x2)6{x3) (see the discussion below eq. (HH)), whose associated 
Fourier transform is g sm{^Y^) / (y) ■ Then, in the exponent of the second 
member in eq. (jl7j) . a typical relevant term giving the leading order behav- 
ior in L is, 



{gmg^ f d^k 1 



27r3 A;2(^^ + a/k) 



sin(fciL/2) 
ki/2 



^ / dkr [kx + a/gl ^ ^ 



where kx = (^2 + ^3)^^^- We see that for small kr (large transverse distances) 
this result is finite, that is, no cut-off of the order of the loop transverse 
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dimension is needed. Therefore, when the complete loop is considered, the 
perimeter law is obtained. This corresponds to string breaking due to the 
creation of light dynamical quark-antiquark pairs out of the vacuum. 

Outside this limiting situations, the quadratic expansion of F is not in 
general a controlled procedure to study infrared behavior. Then, for finite 
masses, it would be interesting to include the effect of nonquadratic terms, 
in order to gain additional information about the scale for the onset of the 
Wilson loop perimeter law. 

In this work, considering the London limit of the SU(2) Georgi-Glashow 
model, we have followed Polyakov's monopole plasma approach, together 
with dual methods for the fermionic matter, to derive in a straightforward 
manner a useful dual representation for the Wilson loop. In this way, a Z(2) 
local symmetry is directly evidenced. This gauging is realized in terms of 
the vector field bosonizing the fermion current, along the U(l) subgroup of 
SU(2) which is unbroken before compactification. 

Although we have considered fermionic quarks, represented by Dirac 
fields, our procedure can be also extended to scalar matter, by following 
the generalized dual methods discussed in ref. pUj . 

The gauging of the Z(N) symmetry has been previously proposed in [TU] , 
using the vortex condensation approach to the SU(N) Georgi-Glashow model 
with scalar matter in the fundamental representation. In that reference, 
a Wilson loop perimeter law finds an heuristic interpretation in terms of 
overlaps. Namely, when Z(N) is local, the Wilson loop operator changes 
the vacuum inside the loop to a gauge equivalent state, while outside the 
loop it is left unchanged. Then, the overlap associated with the operator 
vacuum expectation value is expected to receive nontrivial contributions from 
a boundary effect. 

Here, in the phase where dynamical quarks are coupled, we have seen 
that the theory displays a local Z(2) symmetry, and the relevant source is 
indeed concentrated on a boundary. In this phase, when the large mass limit 
is considered, the scale of localization becomes infinite and the perimeter law 
cannot be established for any given finite loop, as expected. On the other 
hand, in the limit m — 0, an infrared fixed point theory can be used to 
study the Wilson loop behavior. In this perimeter law is obtained, 

corresponding to the creation of light dynamical quark-antiquark pairs. 

Summarizing, string breaking is a widely observed phenomenon in lat- 
tice simulations of (2 + 1)D models. Therefore, it is desirable to obtain a 
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general theoretical understanding of the subject. Here, we have seen that 
the consideration of dual methods to treat the dynamical matter constitutes 
an interesting possibility to deal with compactification and discuss string 
breaking phases coming from the quark sector. 

1 Acknowledgements 

The author would like to thank C. Schat for useful comments. The Conselho 
Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq-Brazil) and the 
Fundagao de Amparo a Pesquisa do Estado do Rio de Janeiro (FAPERJ) are 
acknowledged for the financial support. 

References 

[1] A.M. Polyakov, Phys. Lett. B59, 82 (1975); 
A. M. Polyakov, Nucl. Phys. B 120 (1977) 429. 

[2] G. 't Hooft, Nucl. Phys. B 138 (1978) 1; 

G. 't Hooft, Nucl. Phys. B 190 (1981) 455. 

[3] I. I. Kogan and A Kovner, piep^h/0205026| 

[4] M. N. Chernodub, F. V. Gubarev, M. I. Polikarpov, V. I. Zakharov, 
hep-lat/0 103033; 

N. Seiberg and E. Witten, Nucl. Phys. B 426 (1994) 19; 

R. Auzzi, S. Bolognesi, J. Evslin and K. Konishi, Nucl. Phys. B 673 
(2003) 187. 

[5] I. F. Herbut and B. H. Seradjeh, Phys. Rev. Lett. 91 (2003) 171601; 

H. Kleinert, F. S. Nogueira and A. Sudbo, Phys. Rev. Lett. 88 (2002) 
232001. 

[6] N. E. Mavromatos and S. Sarkar, Int. J. Mod. Phys. B 17 (2003) 251. 

[7] C. DeTar, O. Kaczmarek, F. Karsch and E. Laermann, Phys. Rev. D 59 
(1999) 031501. 



10 



[8] M. N. Chernodub, E. M. Ilgenfritz and A. Schiller, Phys. Lett. B 547 
(2002) 269. 

[9] H. D. Trottier, Phys. Rev. D 60 (1999) 034506. 

[10] C. D. Fosco and A. Kovner, Phys. Rev. D 63 (2001) 045009. 

[11] E. Fradkin and F. A. Schaposnik, Phys. Lett. B 338 (1994) 253; 
CP. Burgess and F. Quevedo, Nucl. Phys. B 421 (1994) 373; 
F. Schaposnik, Phys. Lett. B 356 (1995) 39. 

[12] E. C. Marino, Phys. Lett. B 263 (1991) 63; 

A. Kovner and P. S. Kurzepa, Int. J. Mod. Phys. A 9 (1994) 4669. 

[13] D. G. Barci, C. D. Fosco and L. E. Oxman, Phys. Lett. B 375 (1996) 
267. 

[14] D. G. Barci, L. E. Oxman and S. P. Sorella, Phys. Rev. D 59 (1999) 
105012; 

D. G. Barci and L. E. Oxman, Nucl. Phys. B 580 (2000) 721. 
[15] L. E. Oxman and S. P. Sorella, Phys. Lett. B 531 (2002) 305. 

[16] N. Bralic, E. Fradkin, V. Manias and F. A. Schaposnik, Nucl. Phys. B 

446 (1995) 144; 

J. G. Le Guillou, E. Moreno, C. Nunez and F. A. Schaposnik, Phys. 
Lett. B 409 (1997) 257; 

R. Banerjee, Phys. Lett. B 358 (1995) 297; 

R. Banerjee, Nucl. Phys. B 465 (1996) 157. 

[17] D. Antonov, Surveys High Energ. Phys. 14 (2000) 265, hep-th/9909209j 

[18] L. M. Krauss and F. Wilczek, Phys. Rev. Lett. 62 (1989) 1221. 

[19] J. Froehlich, R. Goetschmann and P. A. Marchetti, J. Phys. A 28 (1995) 
1169. 

[20] G. D. Fosco, V. E. R. hemes, L. E. Oxman and S. P. Sorella, O. S. 
Ventura, Annals of Phys. 290 (2001) 27. 



11 



